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ABSTRACT 

Measurements  of  surface  forward  scattering  loss  between  20  and  50  kHz  were  made 
from  the  research  platform  FLIP  in  January  1992.  These  measurements  were  analyzed 
for  residual  loss  (defined  as  the  loss  in  excess  of  spreading  and  chemical  absorption), 
which  is  attributed  to  bubbles  and  referred  to  as  surface  bubble  loss  (SBL).  With 
these  data,  a  new,  semiempirical  model  for  SBL  was  developed  that  improves  upon  a 
model  currently  used.  SBL  is  a  key  input  parameter  for  predicting  bubble  attenuation 
in  simulations  of  weapon  system  performance  in  near-surface  environments. 
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EXECUTIVE  SUMMARY 

The  Applied  Physics  Laboratory  at  the  University  of  Washington  (APL-UW) 
conducted  measurements  of  forward  scattering  loss  between  20  and  50  kHz  from 
the  research  platform  FLIP  off  the  coast  of  California  in  January  1992.  This  work 
was  sponsored  by  the  Office  of  Naval  Technology  with  technical  management  by  the 
Applied  Research  Laboratory,  Pennsylvania  State  University.  The  main  objective  was 
to  obtain  coordinated  measurements  of  acoustic  surface  scattering,  the  subsurface 
bubble  field,  and  the  sea  state  in  order  to 

•  understand  more  fully  the  conditions  imposed  by  the  near-surface  environment 
that  affect  the  performance  of  sonar  signal-processing  algorithms 

•  improve  high-frequency  acoustic  models  used  in  sonar  system  simulation,  mod¬ 
eling,  and  analysis. 

The  surface  forward  scattering  loss  measurements  were  analyzed  for  residual  loss 
(defined  as  loss  in  excess  of  spreading  and  chemical  absorption),  which  is  attributed 
to  bubbles  and  referred  to  as  surface  bubble  loss  (SBL).  A  new,  semiempirical  model 
for  SBL  is  developed  from  the  data  that  improves  upon  a  current  model  for  SBL. 
The  new  model  predicts  the  mean  energy  loss  in  surface  bounce  paths  that  is  due  to 
extinction  from  near-surface  bubbles.  SBL  is  a  key  input  parameter  for  predicting 
bubble  attenuation  for  simulation  and  analysis  of  sonar  system  performance  in  near¬ 
surface  environments. 

The  new  model  is 

SBL(dB)  =  t/>4m/s 

Sin  6/ 

SBL  (dB)  =  SBL1„^,  eio(t^-4)^  U  <  4  m/s  , 

where  U  is  wind  speed  measured  10  m  above  the  sea  surface,  F  is  acoustic  frequency 
in  kilohertz,  and  9  is  the  nominal  grazing  angle  of  the  surface  bounce  path  (9  >  0®). 
The  model  includes  a  wind  speed  threshold  of  4  m/s  to  account  for  the  threshold  of 
breaking  waves  (Beaufort  velocity)  and  subsequent  production  of  bubbles. 
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1.  INTRODUCTION 

This  report  presents  results  of  surface  forward  scattering  loss  measurements  made 
in  the  open  ocean  between  20  and  50  kHz.  It  also  describes  a  new  model  for  predicting 
surface  loss  due  to  bubbles,  based  on  these  measurements. 

The  measurements  were  made  from  the  research  platform  FLIP  off  the  coast  of 
California  in  January  1992.  The  objective  of  the  FLIP  experiment  was  to  obtain 
coordinated  measurements  of  acoustic  surface  scattering,  the  subsurface  bubble  field, 
and  sea  state.  The  measurements  will  be  used  to  (1)  understand  more  fully  the 
conditions  imposed  by  the  near-surface  environment  that  affect  the  performance  of 
sonar  signal  processing  algorithms,  and  (2)  improve  high-frequency  acoustic  models 
used  in  sonar  system  modeling  and  analysis. 

Besides  the  scattering  loss  measurements,  the  experiment  also  included  measure¬ 
ments  of  the  strength  and  coherence  of  monostatic  backscatter  and  bistatic  for¬ 
ward  scatter,  which  will  be  discussed  in  separate  reports.  Some  of  the  data  have 
been  transferred  to  the  Applied  Research  Laboratory  at  Pennsylvania  State  Univer¬ 
sity  (ARL/PSU)  for  processing  into  estimates  of  channel  scattering  function  [1]  and 
bistatic  forward  scatter  [2]. 

The  surface  forward  scattering  loss  measurements  were  analyzed  for  residual  loss 
(defined  as  loss  in  excess  of  spreading  and  chemical  absorption)  attributed  to  resonant 
extinction  from  near  surface  bubbles.  This  will  be  referred  to  as  surface  bubble 
loss  (SBL).  The  SBL  estimates  from  the  FLIP  experiment  were  used  to  develop  a 
semiempirical  predictive  model  for  SBL,  which  predicts  mean  energy  loss  in  surface 
bounce  paths  due  to  extinction  from  near-surface  bubbles.  We  recommend  that  this 
new  model  replace  the  SBL  model  currently  used.  SBL  is  a  key  input  parameter 
for  predicting  bubble  attenuation  in  simulations  of  weapon  systems  in  near-surface 
environments. 

This  work  is  part  of  the  Office  of  Naval  Technology’s  Exploratory  Development 
Program  with  technical  management  by  ARL/PSU.  The  ONT  (now  ONR-T)  program 
is  primarily  concerned  with  exploratory  development  within  the  frequency  range  of  10 
to  500  kHz,  with  applications  directed  toward  sonar  signal  processor  and  performance 
analysis  technologies.  A  component  of  the  FLIP  experiment  was  also  carried  out  as 
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part  of  The  Technical  Cooperation  Program  (TTCP).  Comparative  measurements 
of  the  ocean  bubble  layer  were  made  by  APL-UW,  representing  the  U.S.,  and  the 
Institute  of  Ocean  Sciences  (lOS),  representing  Canada.  lOS  measured  the  bubble 
field  using  its  SEASAT  instrument  package  deployed  from  USNS  De  Steiguer  (T- 
AGOR-12)  which  kept  station  nearby  during  9  days  of  the  18-day  experiment.  Joint 
analysis  of  the  APL  and  lOS  measurements  is  continuing;  preliminary  results  of  this 
work  were  presented  last  fall  [3]. 

Section  2  describes  the  FLIP  experimental  measurements.  Section  3  reviews  the 
method  used  for  estimating  surface  bubble  loss,  and  Section  4  summarizes  the  mod¬ 
eling  methods.  A  discussion  and  summary  are  presented  in  Section  5.  Measurement 
errors  are  summarized  in  Appendix  A.  Additional  model/data  comparisons  are  con¬ 
tained  in  Appendix  B. 
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2,  EXPERIMENTAL  DESCRIPTION 

Figure  1  shows  the  track  of  FLIP  during  the  experiment.  The  measurements 
discussed  in  this  report  were  made  between  18  and  24  January.  During  this  period, 
wind  speeds  ranged  from  approximately  1  to  9  m/s,  and  the  sea  often  had  a  large 
swell  component  with  a  period  of  ~11  s  and  significant  wave  height  {Hi/ 3)  between 
1  and  2  m. 

The  geometry  used  for  the  surface  forward  scattering  measurements  is  shown  in 
Fig.  2.  Signals  were  transmitted  from  one  of  three  ITC-1032  transducers  suspended 
from  the  spar  buoy  at  depths  of  27,  57,  and  147  m.^  The  spar  buoy  was  tethered  to 
FLIP  by  a  1000-m  cable  and  ranged  from  500  to  1000  m  from  the  vessel. 

CW  pulses  (1,  4,  or  8  ms)  were  transmitted  at  1-s  intervals  in  a  four-pulse  sequence, 
during  which  the  frequency  changed  from  20  to  50  kHz  in  10-kHz  steps.  A  typical 
run  consisted  of  200  pings,  equivalent  to  50  pings  of  a  single  frequency  transmitted 
at  4-s  intervals. 

The  signals  were  received  at  FLIP  using  a  second  ITC-1032  transducer  located  on 
the  tip  of  the  13-m  subsurface  boom  attached  to  FLIP’S  hull  at  a  depth  of  28.5  m. 
The  boom  was  designed  for  vertical-incidence  measurements  of  the  bubble  layer  but 
also  made  an  excellent  receiver  location,  since  surface  scatter  could  be  distinguished 
from  rescatter  from  the  FLIP’S  hull,  which  arrived  approximately  17  ms  after  the 
direct  arrival. 

Environmental  measurements  consisted  of  mean  wind  speed  and  direction  (10-min 
averages)  taken  10  and  25  m  above  the  sea  surface,  surface  wave-height  spectra,  air 
and  sea  temperature,  and  conductivity-temperature-depth  (CTD)  casts  (twice  daily, 
typically  at  0800  and  2000).  Additional  details  on  the  environmental  measurements 
are  available  in  Ref.  4. 


^ These  depths  varied  ±0.5  m  owing  to  heave  motion  of  the  buoy. 


Figure  1.  FLIP  drift  track  as  determined  from  the  GPS  navigation  system.  Loca¬ 
tions  at  0000  hours  local  time  are  indicated  along  the  track. 
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Figure  2.  Experimental  geometry  for  surface  forward  scattering  loss  measurements. 
The  spar  buoy  range  varied  between  500  and  1000  m. 
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3.  DATA  ANALYSIS  METHODS 
3.1  Preliminary  Discussion 

In  general  a  signal  forward  scattered  from  the  sea  surface  experiences  both  time 
spreading  (pulse  elongation)  and  energy  loss.  The  latter  is  attributed  to  the  extin¬ 
guishing  effects  of  the  near-surface  bubble  field,  through  which  surface  bounce  paths 
must  traverse.  This  effect  is  known  as  surface  bubble  loss  (SBL).  Reference  5  gives 
a  model  for  the  time  dependence  of  forward-scattered  intensity.  This  model  pre¬ 
dicts  the  build-up  and  decay  of  mean  intensity  as  a  function  of  geometry  and  surface 
roughness  (assuming  high  frequency  linoit,  Gaussian  surface  statistics).  The  model 
matches  closely  the  ensemble  average  of  simulated  forward-scattered  pulses  using 
specular  point  theory  [6]  as  well  as  data  from  FLIP. 

In  the  absence  of  attenuation  from  near-surface  bubbles,  the  time  integrated  in¬ 
tensity  can  be  modeled  with  a  0  dB  loss,  this  result  having  been  explicity  proven 
using  the  Kirchhoff  approximation  and  confirmed  with  experimental  data  [3,7].  We 
thus  estimate  SBL  by  estimating  the  residual  loss  (in  decibels)  in  total  energy,  after 
accounting  for  spreading  loss  and  chemical  absorption,  the  latter  being  temperature 
and  salinity  dependent.  The  SBL  estimate  is  equivalent  to  the  residual  loss  estimate. 

A  necessary  condition  for  the  0  dB  loss  is  that  the  transmit  and  receive  beamwidths 
be  greater  than  the  angular  width  of  the  surface-scattered  intensity.  In  practice,  this 
criterion  is  met  with  vertical  beamwidths  of  20-25°  and  horizontal  beam  widths  of 
5-10°  [5].  For  the  FL/P  experiment,  the  transmit  and  receive  beams  were  omnidirec¬ 
tional  to  avoid  any  effects  caused  by  narrow  beams. 

3.2  SBL  Estimation  Procedure 

As  noted  above,  pulse  energy  is  generally  spread  in  time.  In  the  absence  of  a 
bubble  layer,  however,  the  total  energy  level  scattered  from  the  surface,  Eq,  can  be 
computed  by 

Eo  (dB)  =  5L  -  TL  +  10  log  r  ,  (1) 

where  TL  is  transmission  loss,  which  includes  geometric  spreading  and  chemical  ab¬ 
sorption,  SL  is  peak  source  level,  and  r  is  the  pulse  length. 
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For  each  run,  a  ray  analysis  is  performed  using  the  nearest-time  CTD  cast  for  the 
sound  speed  profile.  From  the  ray  analysis,  we  get  direct  and  surface-path  eigenrays 
that  match  the  measured  travel  times,  surface  path-spreading  loss  (typically  20  log 
of  the  range),  and  nominal  surface  grazing  angle  (6).  Chemical  absorption  losses  are 
estimated  from  [8],  based  on  the  temperature  and  salinity  at  30  m. 

The  output  of  our  acoustic  receiving  system  is  a  signal  mixed  down  to  5  kHz,  with 
a  nominal  bandwidth  of  10  kHz.  The  5-kHz  signal  is  digitized  at  a  20-kHz  rate,  and 
later  digitally  bandpass  filtered  between  4  and  6  kHz  to  remove  distortions  and  noise 
that  may  have  been  introduced  during  the  conditioning  phase.  Echo  envelopes  are 
computed  from  the  Hilbert  transform  of  each  ping.  We  form  an  equivalent  estimate 
of  energy  from  the  data,  Ei,  by  integrating  the  intensity  over  the  time  span  of  each 
surface  arrival  (see  Figs.  3  and  4).  Note  that  the  leading  edges  of  the  direct  path  are 
first  aligned  to  remove  effects  due  to  source  and  receiver  motion.  The  integration  is 
truncated  approximately  17  ms  after  the  start  of  direct  arrival  to  avoid  integrating 
rescatter  from  FLIP'S  hull.^  For  each  ping  an  estimate  of  the  residual  loss  is  computed 
which  is  equal  to  Eq  —  Ei.  The  residual  loss  is  our  estimate  of  the  surface  bubble 
loss.  The  mean  loss  (SBL)  is  computed  from  the  energy  ensemble  average  for  each 
run;  i.e.,  (SBL)  =  Eq  —  (F^i),  where  [E-^)  represents  time  integration  of  the  ensemble 
average. 

In  runs  for  which  the  direct  and  surface  paths  overlap,  it  is  not  possible  to  integrate 
the  surface  path  separately.  Here,  we  use  an  alternate  method  [9]  for  estimating  the 
surface  bubble  loss  based  on  fitting  the  ensemble  mean  of  surface-scattered  intensity 
using  the  pulse  elongation  model  in  [5].  The  pulse  elongation  model  requires  an 
estimate  of  rms  slope,  and  the  surface  loss  is  then  adjusted  to  fit  the  data. 

The  two  SBL  estimation  procedures  are  illustrated  in  Figs.  3  and  4  using  data 
from  run  114  made  during  conditions  characterized  by  fairly  frequent  white  capping 
(Beaufort  Wind  Scale  4-5)  with  U  =  7  m/s.  To  fit  the  data  ensemble  mean  with 
the  pulse  elongation  model,  the  30-kHz  data  (Fig.  3)  requires  3  dB  of  SBL,  while 


^Typically  less  than  1%  of  the  energy  arrives  after  this  point,  amounting  to  negligible  error. 
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time,  ms 


Figure  3.  Surface  forward- scattered  intensity  vs  time  for  30-kHz  run  114,  showing 
ensemble  mean  (*)  and  ping  number  10  (o).  The  dashed  line  represents  the  time 
spread  model  in  Ref.  5  fitted  to  the  ensemble  data;  in  this  case,  the  model  is  given 
an  SBL  input  of  3  dB  in  order  to  fit  the  data. 

the  40-kHz  data  (Fig.  4)  requires  4.5  dB  of  SBL.  These  estimates  match  closely  those 
derived  by  direct  integration  (3.0  and  4.6  dB).  In  each  figure,  ping  number  10  for  each 
frequency^  is  shown  including  the  region  of  integration  for  computing  Ei,  starting  at 
the  beginning  of  the  surface  arrival  and  ending  just  before  rescatter  from  FLIP'S  hull 
commences  (~17  ms). 


^Each  30  kHz  ping  is  transmitted  1  s  before  each  40  kHz  ping. 
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Figure  4.  Surface  forward  scattered  intensity  vs  time  for  40-kHz  run  114  showing 
ensemble  mean  (*)  and  ping  number  10  (o).  The  dashed  line  represents  the  time 
spread  model  in  Ref,  5  fitted  to  the  ensemble  data;  in  this  case,  the  model  is  given 
an  SBL  input  of  4.5  dB  in  order  to  fit  the  data. 


3.3  Results 

Table  I  summarizes  the  SBL  database  derived  from  the  FLIP  experiment.  When 
possible,  the  direct  integration  method  has  been  used  to  estimate  SBL  since  with 
this  method  we  can  also  estimate  ping-to-ping  fluctuations  in  SBL,  For  example,  the 
30-kHz  ping  number  10  experienced  a  6.5  dB  SBL  compared  with  an  ensemble  mean 
loss  of  3  dB,  and  the  40-kHz  ping  number  10  experienced  a  10  dB  loss  compared  with 
an  ensemble  mean  loss  of  4.5  dB.  Clearly,  fluctuations  in  SBL  can  be  large.  For  a 
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Table  I.  Summary  of  SBL  database  derived  from  FLIP  experiment. 


EXP. 

RUN 

RANGE 

(m) 

FREQ. 

(kHz) 

6 

(deg) 

WIND"^ 

(m/s) 

(SBL) 

(dB) 

1 

109 

874 

20 

10.9 

8.0 

1.2 

2 

109 

874 

30 

10.9 

8.0 

3.7 

3 

110 

877 

20 

10.9 

7.8 

1.0 

4 

no 

874 

30 

10.9 

7.8 

3.8 

5 

no 

874 

40 

10.9 

7.8 

5.2 

6 

no 

874 

50 

10.9 

7.8 

3.5 

7 

111 

855 

20 

4.5 

7.8 

2.7 

8 

111 

855 

30 

4.5 

7.8 

6.3 

9 

111 

855 

40 

4.5 

7.8 

9.6 

10 

111 

855 

50 

4.5 

7.8 

9.0 

11 

114 

995 

20 

11.1 

7.0 

1.2 

12 

114 

995 

30 

11.1 

7.0 

3.0 

13 

114 

995 

40 

11.1 

7.0 

4.6 

14 

114 

995 

50 

11.1 

7.0 

4.0 

15 

115 

995 

20 

11.1 

5.3 

1.2 

16 

115 

995 

30 

11.1 

5.3 

2.4 

17 

115 

995 

40 

11.1 

5.3 

2.6 

18 

115 

995 

50 

11.1 

5.3 

1.6 

19 

116 

540 

20 

17.6 

4.1 

1.7 

20 

116 

540 

30 

17.6 

4.1 

0.2 

21 

116 

540 

40 

17.6 

4.1 

0.1 

22 

116 

540 

50 

17.6 

4.1 

0.2 

23 

117 

540 

20 

17.6 

4.1 

1.7 

24 

117 

540 

30 

17.6 

4.1 

0.8 

25 

117 

540 

40 

17.6 

4.1 

0.4 

26 

117 

540 

50 

17.6 

4.1 

0.5 

27 

118 

532 

20 

9.1 

4.1 

0.5 

28 

118 

532 

30 

9.1 

4.1 

0.5 

29 

118 

532 

40 

9.1 

4.1 

2.6 

30 

118 

532 

50 

9.1 

4.1 

3.0 

31 

132 

1000 

20 

9.3 

6.2 

-0.8 

32 

132 

1000 

30 

9.3 

6.2 

1.7 

33 

132 

1000 

40 

9.3 

6.2 

1.5 

34 

132 

1000 

50 

9.3 

6.2 

2.5 

35 

133 

668 

20 

14.3 

4.3 

-0.1 

36 

133 

668 

30 

14.3 

4.3 

0.2 

37 

133 

668 

40 

14.3 

4.3 

-0.2 

38 

133 

668 

50 

14.3 

4.3 

0.2 

39 

138 

666 

20 

14.5 

4.4 

0.7 

40 

138 

666 

30 

14.5 

4.4 

1.3 

41 

138 

666 

40 

14.5 

4.4 

2.5 

42 

138 

666 

50 

14.5 

4.4 

0.2 

43 

139 

666 

20 

14.5 

4.8 

0.7 

44 

139 

666 

30 

14.5 

4.8 

1.3 

45 

139 

666 

40 

14.5 

4.8 

1.5 

46 

139 

666 

50 

14.5 

4.8 

0.2 

‘^Measured  10  m  above  sea  surface. 
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useful  measure  of  the  size  of  these  fluctuations  we  adapt  an  expression  for  transmission 
loss  studies  [10]  to  our  estimates  of  surface  loss  due  to  bubbles  (SBL)  and  define 

SBL+  =  mog[Eo/{{E^)-a)] 


and 

SBL_  =  101og[Eo/((^i)  +  (7)]  , 

where  a  is  the  standard  deviation  of  the  time- integrated  intensity;  SBL+  represents  a 
typical  high  loss  value,  and  SBL_  represents  a  low  loss  value.  We  expect  the  majority 
of  ping-to-ping  fluctuations  in  surface  loss  to  lie  between  these  bounds.  From  the 
FLIP  data,  the  low  and  high  loss  bounds  are  such  that  SBL_  w  (SBL)  —  3  dB  and 
SBL+  Fs  (SBL)  -f  5  dB.  The  data  suggest  that  fluctuations  in  loss  estimates  tend  to 
be  dominated  by  fluctuations  associated  with  incoherent  scattering  from  a  random 
surface  [11].  Under  these  conditions,  it  is  very  difficult  to  resolve  the  component  of 
variation  due  solely  to  the  bubble  field. 
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4.  MODEL  FOR  SURFACE  BUBBLE  LOSS 

4.1  Preliminary  Discussion 

The  two-way  transmission  loss  attributable  to  bubbles  [12]  can  be  generalized  to 
an  arbitrary  grazing  angle,  giving  the  following  expression  for  SBL  [5]. 

SBL  (dB)  =  f  s,{z)dz  ,  (2) 

smt^  JD 

where  6  is  the  nominal  grazing  angle  for  surface  arrivals;  Se{z)  is  the  total  extinction 
cross  section  per  unit  volume,  which  is  integrated  over  the  bubble  layer  thickness  D. 
The  bubble  layer  thickness  in  this  case  corresponds  to  the  average  entrainment  depth 
of  resonant  bubbles. 

We  use  Eq.  (2)  as  a  guide  for  developing  a  predictive  model  for  SBL  as  a  func¬ 
tion  of  grazing  angle  6  and  frequency,  with  wind  speed  as  an  environmental  predictor 
of  both  bubble  concentration  and  entrainment  depth.  The  l/sin0  dependence  is 
well  supported  by  the  FLIP  database  and  other  measurements  [13],  but  modeling 
the  depth-integrated  Se(2)  is  more  difficult.  In  principle,  one  can  infer  the  necessary 
depth-integrated  properties  of  SeC-^)  from  vertical-incidence  sonar  measurements.  But 
predicting  bubble  attenuation  effects  from  bubble  scattering  data  have,  to  date,  given 
inconsistent  results  [13,14].  The  source  of  error  may  be  in  the  bubble  damping  con¬ 
stant,  or  perhaps  in  assumptions  regarding  the  horizontal  and  vertical  variations  in 
bubble  density. 

4.2  Modeling  Surface  Bubble  Loss 

The  approach  we  take  is  to  use  the  actual  surface  loss  data  to  estimate  parameters 
in  a  semiempirical  model  of  the  form 

SBL  (dB)  =  ,  !/  >  4  m/s  (3) 

sm  6^ 

SBL  (dB)  =  SBLlu^,  e^o([7-4)  ^  t/  <  4  m/s  , 

where  U  is  wind  speed  measured  10  m  above  the  sea  surface,  F  is  acoustic  frequency 
in  kilohertz,  and  parameters  pi  to  ps  are  derived  from  the  data.  We  assume  bubble 
production  and  entrainment  follows  a  wind  speed  power  law,  consistent  with  well 
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established  power  law  forms  for  whitecap  coverage  [15];  from  typical  bubble  size 
spectra  [16],  we  would  also  expect  the  frequency  dependence  to  be  a  power  law.  We 
retain  the  1/  sin  0  dependence  to  be  consistent  with  Eq.  (2),  and  include  a  wind  speed 
threshold  of  4  m/ s  to  account  for  the  breaking  wave  threshold,  or  Beaufort  velocity, 
often  reported  to  be  between  2  and  5  m/s.  For  the  FLIP  measurements,  we  observed 
essentially  no  surface  loss  for  wind  speeds  less  than  4  m/s.  The  transition  to  0  dB 
loss  is  approximated,  primarily  for  continuity  in  simulations,  by  a  steep  exponential 
decay  such  that  the  loss  at  3.5  m/s  is  about  1%  of  the  loss  at  4  m/s  (i.e.,  SBL|u^^). 

The  three  parameters  are  estimated  by  fitting  a  nonlinear  curve  to  Eq.  (3),  using 
the  FLIP  data  in  Table  I  (freq,  U,  0,  and  (SBL)).  The  Nelder-Meade  algorithm^  is 
used  to  minimize  the  total  absolute  value  error  between  model  and  data  vectors.  The 
fitted  model  parameters  are 

Pi  =  1.26  ±  0.48  X  10-® 

P2  =  1.57  ±0.23 
P3  =  0.85  ±  0.22  , 

where  the  ±  values  represent  the  95%  confidence  interval  for  the  parameter  estimates, 
as  generated  by  the  bootstrap  algorithm  [17]. 


“^See  MATLAB  User’s  Guide. 
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5.  DISCUSSION  AND  SUMMARY 

Figures  5-8  show  representative  model/data  comparisons  between  20  and  50  kHz. 
To  show  data  collected  at  different  grazing  angles,  all  data  have  been  normalized  to 
a  10°  grazing  angle  for  this  comparison,  using  the  formula 

SBL(10“)  =  4^SBL(«)  ,  (4) 

sm  10° 

where  SBL(0)  is  the  actual  measured  loss®  at  grazing  angle  $. 

The  scatter  in  the  data  generally  increases  for  smaller  wind  speeds,  owing  in  part  to 
measurement  uncertainties  (approximately  ±2.5  dB;  see  Appendix  A).  Nevertheless, 
the  model  provides  a  reasonable  representation  of  the  data,  with  nearly  all  data  points 
falling  within  ±2  dB  of  the  predicted  curve.  The  sharp  drop  in  each  predicted  curve 
at  a  wind  speed  equal  to  4  m/s  represents  the  wind  speed  threshold  for  the  onset 
of  breaking  waves  and  subsequent  production  of  bubbles.  The  entire  model  is  now 
summarized  with  parameters 

SBL  =  t/>4m/s  (5) 

sin  6/ 

SBL  =  SBL|„^,  eio(i/-4)^  (7  <  4  m/s  , 

where  SBL  is  in  decibels.  Note  that  grazing  angle  $  must  be  restricted  to  ^  >  0°. 

Using  the  known  environmental  conditions  and  fixed  experimental  geometries  from 
these  two  experiments,  the  model  describes  the  FLIP  measurements  to  within  ±2  dB, 
with  a  total  rms  error  of  approximately  1  dB.  Fluctuations  in  model  estimates  using 
the  bootstrap  algorithm  gave  similar  values.  In  general,  however,  uncertainties  in 
environmental  conditions  and  run  geometry,  e.g.,  the  grazing  angle  of  the  surface 
bounce  path,  will  tend  to  magnify  errors  in  model  prediction,  and  we  set  a  more 
realistic  bound  for  model  uncertainty  to  ±3  dB.  Additional  model/data  comparisons 
using  SBL  measurements  from  two  other  experiments  are  shown  in  Appendix  B. 

The  SBL  model  predicts  the  mean  energy  loss  in  surface  bounce  paths  due  to 
extinction  from  near-surface  bubbles,  where  the  energy,  or  time-integrated  intensity, 


®The  three  negative  values  in  Table  I  have  been  set  to  zero  for  this  comparison. 
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Figure  5.  Comparison  of  SBL  model  with  measured  20-kHz  data  from  FLIP.  Data 
from  various  grazing  angles  have  been  normalized  to  10°  grazing  angle  according  to 
Eq.  (5).  Upper  and  lower  lines  represent  ±2  dB  interval  about  the  predicted  curve. 


Figure  6,  Comparison  of  SBL  model  with  measured  30-kHz  data  from  FLIP.  Data 
from  various  grazing  angles  have  been  normalized  to  10°  grazing  angle  according  to 
Eq.  (5).  Upper  and  lower  lines  represent  ±2  dB  interval  about  the  predicted  curve. 
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Figure  7.  Comparison  of  SBL  model  with  measured  40-kHz  data  from  FLIP.  Data 
from  various  grazing  angles  have  been  normalized  to  10°  grazing  angle  according  to 
Eq.  (5).  Upper  and  lower  lines  represent  ±2  dB  interval  about  the  predicted  curve. 


Figure  8.  Comparison  of  SBL  model  with  measured  50-kHz  data  from  FLIP.  Data 
from  various  grazing  angles  have  been  normalized  to  10°  grazing  angle  according  to 
Eq.  (5).  Upper  and  lower  lines  represent  ±2  dB  interval  about  the  predicted  curve. 
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can  be  modeled  with  a  0  dB  loss  in  the  absence  of  bubbles.  Model  performance  will 
be  degraded  in  conditions  under  which  the  0  dB  loss  case  cannot  be  assumed.  Some 
examples  where  this  could  occur  are  as  follows: 

•  Model  used  to  simulate  data  gathered  with  narrow  transmit  and  receive  beam 
patterns  (see  Section  3.1). 

Data  losses  will  be  greater  than  simulated  SBL  because  of  an  undersampling  of 
the  angular  distribution  of  scattered  intensity  in  the  data. 

•  SBL  model  used  to  simulate  loss  in  the  peak  value  of  short-pulse  data. 

A  short  pulse  is  defined  as  one  with  pulse  duration  r  much  less  than  the  char¬ 
acteristic  pulse  elongation  time  T  for  surface-scattered  arrivals.  In  such  cases, 
the  intensity  will  not  have  reached  steady  state,  and  the  total  energy  will  be 
underestimated  when  based  on  estimates  of  peak  intensity,  resulting  in  data 
losses  greater  than  the  simulated  losses.  Note  that  ideally  the  time  interval 
of  intensity  integration  should  be  2-3  times  the  characteristic  pulse  elongation 
time  T  [5].  A  rough  guide  is  T  =  4.6i?,  where  T  is  in  milliseconds  and  range  R 
is  in  kilometers. 

•  SBL  model  used  to  simulate  conditions  characterized  by  a  highly  directional 
surface  wave  spectrum.  In  this  case,  the  distribution  of  surface  slopes  can  no 
longer  be  considered  isotropic,  and  the  model  error  will  depend  on  the  exact 
orientation  of  the  surface  bounce  paths  with  respect  to  the  surface  waves. 

Expected  uncertainties  in  the  new  SBL  model  are  summarized  as  follows: 

•  The  uncertainty  in  predicted  mean  SBL  is  ±3  dB. 

•  The  ping-to-ping  fluctuations  in  SBL  will  be  approximately  bounded  by  —3 
and  -|-5  dB  about  the  predicted  mean  SBL  value,  based  on  the  SBL+  and  SBL_ 
estimates  discussed  in  Section  3.3. 


TR9307  17 


UNIVERSITY  OF  WASHINGTON*  APPLIED  PHYSICS  LABORATORY, 


REFERENCES 

[1]  Johnson,  B.  L.,  “Preliminary  scattering  function  estimates  from  the  FL/P  experi¬ 

ment,”  TM92-173,  Applied  Research  Laboratory,  Pennsylvania  State  University, 
September  1992. 

[2]  Lang,  D.  C.,  and  R.  L.  Culver,  “A  High  Frequency  Bistatic  Ocean  Surface  Scat¬ 

tering  Strength,”  TM92-342,  Applied  Research  Laboratory,  Pennsylvania  State 
University,  December  1992. 

[3]  Dahl,  P.  H.,  M  V.  Trevorrow,  and  A.  T.  Jessup,  “Comparative  acoustic  measure¬ 

ments  of  the  ocean  bubble  layer,”  J.  Acoust.  Soc.  Am.,  92 (i):  Pt.  2,  1248 
(presented  at  the  Fall  1992  meeting  of  the  Acoustical  Society  of  America). 

[4]  Dahl,  P.  H.,  and  A.  T.  Jessup,  “Bubble  Plumes  and  Breaking  Waves:  Measure¬ 

ments  from  the  R/P  FLIP,  January  1992,”  APL-UW  TM2-92  Applied  Physics 
Laboratory,  University  of  Washington,  May  1992. 

[5]  Thorsos,  E.  I.,  “Surface  Forward  Scattering  and  Reflection,”  APL-UW  7-83,  Ap¬ 

plied  Physics  Laboratory,  University  of  Washington,  May  1984. 

[6]  Funk,  D.,  and  K.  Williams,  “A  physically  motivated  simulation  technique  for  high- 

frequency  forward  scattering  derived  using  specular  point  theory,”  J.  Acoust. 
Soc.  Am.,  91,  2606-2614,  1991. 

[7]  McDonald,  J.  F.,  and  R.  C.  Spindel,  “Implications  of  Fresnel  corrections  in  a 

non-Gaussian  surface  scatter  channel,”  J.  Acoust.  Soc.  Am.,  50,  746-757,  1971. 

[8]  Francois,  R.  E.,  and  G.  R.  Garrison,  “Sound  absorption  based  on  ocean  mea¬ 

surements.  Part  II:  Boric  acid  contribution  and  equation  for  total  absorption,” 
J.  Acoust.  Soc.  Am.,  12,  1879-1890,  1982. 

[9]  Thorsos,  E.  I.,  “High  frequency  surface  forward  scattering  measurements,”  pre¬ 

sented  at  the  108th  meeting  of  the  Acoustical  Society  of  America,  October  1984. 

[10]  Chapman,  D.  M.  F.,  “A  simple  estimate  of  propagation  loss  fluctuations  due  to 
modal  interference,”  J.  Acoust.  Soc.  Am.,  85,  1097-1106,  1989. 


18  TR9307 


UNIVERSITY  OF  WASHINGTON*  APPLIED  PHYSICS  LABORATORY. 


[11]  Winebrenner,  D.  P.,  “Acoustic  Backscattering  from  Sea  Ice  at  10-100  kHz,”  APL- 
UW  TR9017,  Applied  Physics  Laboratory,  University  of  Washington,  January 
1991. 

[12]  Medwin,  H.,  “In  situ  acoustic  measurements  of  microbubbles  at  sea,”  J.  Geophys. 
Res.,  82,  971-976,  1977. 

[13]  McConnell,  S.  O.,  and  P.  H.  Dahl,  “Vertical  incidence  backscatter  and  surface 
forward  scattering  from  near-surface  bubbles,”  Proceedings,  Oceans  ’91,  1991. 

[14]  Medwin,  H.,  “In  situ  acoustic  measurements  of  bubble  populations  in  coastal 
waters,”  J.  Geophys.  Res.,  75,  8239-8248,  1970. 

[15]  Monahan,  E.  C.,  “Occurrence  and  evolution  of  acoustically  relevant  sub-surface 
bubble  plumes  and  their  associated,  remotely  monitorable,  surface  whitecaps,” 
Natural  Sources  of  Underwater  Sound,  503-517,  B.  R.  Kerman  (ed.),  Kluwer 
Academic,  the  Netherlands,  1993. 

[16]  Medwin,  H.,  and  N.  Breitz,  “Ambient  and  transient  bubble  spectral  densities  in 
quiescent  seas  and  under  spilling  breakers,”  J.  Geophys.  Res.,  94,  12,751-12,759, 
1989. 

[17]  Efron,  B.,  and  R.  Tibshirani,  “Bootstrap  methods  for  standard  errors,  confidence 
intervals,  and  other  measures  of  statistical  accuracy,”  Statistical  Science,!,  54-77, 
1986. 

[18]  McConnell,  S.  O.,  “Acoustic  measurements  of  bubble  densities  at  15-50  kHz,”  Sea 
Surface  Sound:  Natural  Mechanisms  of  Surface  Generated  Noise  in  the  Ocean, 
237-252,  B.  R.  Kerman  (ed.),  Kluwer  Academic,  the  Netherlands,  1988. 


TR9307  19 


UNIVERSITY  OF  WASHINGTON*  APPLIED  PHYSICS  LABORATORY, 


APPENDIX  A 

Calibration  and  Measurement  Errors 

The  transmitting  system  (three  ITC  projectors  and  1000  m  of  EM  cable)  and  re¬ 
ceiving  systems  were  calibrated  at  the  APL-UW  calibration  facility  during  the  month 
before  the  experiment;  overall  accuracy  of  the  calibration  system  is  ±1.5  dB. 

In  the  course  of  analyzing  the  data,  however,  we  observed  a  systematic  bias, 
particularly  for  runs  made  at  50  kHz.  For  example,  measurements  gathered  under 
conditions  characterized  by  wind  speeds  ~4  m/s  and  the  absence  of  wavebreaking 
activity  for  several  hours  are  expected  to  show  a  0  dB  loss.  This  was  the  case  for  the 
30-kHz  runs  (within  ±2  dB),  but  the  50-kHz  runs  showed  a  systematic  average  nega¬ 
tive  loss  of -4  dB,  while  the  20-kHz  runs  showed  a  small  average  positive  loss  of  1.2  dB. 
We  attribute  this  result  to  differences  in  cable  loading  during  the  calibration  (cable 
spooled)  and  during  the  in  situ  measurements  (cable  laid  out  in  seawater).  It  is  not 
possible  to  reconstruct  the  exact  in  situ  complex  impedance  of  the  cable.  But  simple 
modeling  of  the  cable/transducer  system  using  a  range  of  possible  lumped  circuit  pa¬ 
rameters  reproduced  a  similar  trend  for  the  frequencies  farthest  from  the  transducer 
resonance  (20  and  50  kHz).  Thus,  based  on  five  in  situ,  bubble-free  measurements 
of  surface  loss,  we  have  adjusted  the  source  level  from  what  our  calibrations  would 
normally  give  in  the  following  way:  20  kHz  (down  by  1.2  dB),  30  kHz  (no  change), 
40  kHz  (no  change),  50  kHz  (up  by  4.0  dB).  Similar  estimates  of  the  source  level  bias 
were  made  using  one  run  with  an  extremely  stable  direct  path. 

Another  source  of  error  is  the  TL  estimate  in  Eq.  (1),  which  is  divided  into  (1) 
chemical  absorption  error  (±1.5  dB)  and  (2)  spreading  loss  error  (±1  dB).  The  error 
component  due  to  absorption  is  based  on  the  expected  variation  (±10%)  for  the 
temperature-  and  salinity-dependent  absorption  coefficient,  evaluated  at  50  kHz  and 
using  the  maximum  range  of  our  measurements  (1000  m).  The  error  component  due 
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to  spreading  is  caused  by  differences  in  sound  speed  conditions  between  the  actual 
experimental  run  and  what  was  measured  during  the  CTD  cast.  As  a  guide,  we 
use  the  average  difference  between  expected  spreading  loss  for  isovelocity  conditions 
(approximately  20  log  of  the  range)  and  the  ray  trace  estimate,  setting  this  error  to 
±1  dB. 

The  calibration  errors,  errors  in  the  TL  estimate,  and  errors  due  to  sample  size 
(less  than  1  dB  in  our  data)  sum  independently,  giving  approximately  ±2.5  dB  for 
an  estimate  of  the  total  error  that  applies  independently  to  each  experimental  run. 
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APPENDIX  B 

Additional  Model/Data  Comparisons 

The  SBL  measurements  from  the  FT/P  experiment  were  made  under  open  ocean, 
unlimited  fetch  conditions,  with  the  received  signal  consisting  entirely  of  direct  and 
surface  bounce  paths.  They  represent  the  largest  set  of  combined  acoustic  and  sea 
state  measurements  from  which  to  generate  a  semiempirical  model  for  SBL.  A  smaller 
set  of  SBL  measurements,  taken  in  shallow  water  conditions,  was  obtained  in  the 
Quinault  [9]  and  Whidbey  Island  [13,18]  experiments.  Model  comparisons  using  these 
data  in  addition  to  the  FLIP  data  are  shown  in  Figs.  B1-B3.  Note  that  as  in  Figs. 
5-8  data  have  been  normalized  to  10®  grazing  to  facilitate  comparison. 

In  the  Quinault  experiment,  SBL  was  measured  in  continental  shelf  waters  (depth 
50  m)  approximately  17  km  off  the  Washington  coast.  In  terms  of  sea  surface  proper¬ 
ties,  this  site  can  be  considered  open  ocean.  In  the  Whidbey  Island  experiment,  SBL 
was  measured  in  inland  waters  (depth  30  m)  near  the  junction  of  the  Strait  of  Juan 
de  Fuca  and  Puget  Sound.  Sea  surface  properties  at  this  site  are  strongly  dependent 
upon  wind  direction.  In  both  these  experiments  the  received  signal  consisted  of  a 
mixture  of  surface  and  bottom  paths.  The  Quinault  measurements  showed  evidence 
of  very  high  losses  that  in  some  cases  did  not  correspond  with  wind  speed.  For  the 
Whidbey  Island  data  the  severe  fetch  limitation  may  be  a  factor  in  determining  bub¬ 
ble  dispersion  and  SBL,  and  uncertainties  in  wind  speed  measurements  may  be  a 
source  of  error. 
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Figure  Bl.  Comparison  of  SBL  with  measured  20-kHz  data  from  FLIP  (open  cir¬ 
cles),  Quinault  (stars),  and  Whidbey  Island  (crosses).  Data  from  various  grazing 
angles  have  been  normalized  to  10°  grazing  angle  according  to  Eq.  (5).  Upper  and 
lower  lines  represent  i3  dB  interval  about  the  predicted  curve. 


Figure  B2.  Comparison  of  SBL  with  measured  30-kHz  data  from  FLIP  (open  cir¬ 
cles),  Quinault  (stars),  and  Whidbey  Island  (crosses).  Data  from  various  grazing 
angles  have  been  normalized  to  10°  grazing  angle  according  to  Eq.  (5).  Upper  and 
lower  lines  represent  i3  dB  interval  about  the  predicted  curve. 
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Figure  B3.  Comparison  of  SBL  with  measured  40-kHz  data  from  FLIP  (open  cir¬ 
cles),  and  Whidbey  Island  (crosses).  Data  from  various  grazing  angles  have  been 
normalized  to  10°  grazing  angle  according  to  Eq.  (5).  Upper  and  lower  lines  repre¬ 
sent  ±3  dB  interval  about  the  predicted  curve. 
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